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Qh. Abstract 



Most models designed to study the bidirectional movement of cargos as they are driven 
by molecular motors rely on the idea that motors of different polarities can be coordinated 
If^ • by external agents if arranged into a motor-cargo complex to perform the necessary work [Ij. 

Although these models have provided us with important insights into these phenomena, there 
are still many unanswered questions regarding the mechanisms through which the movement 
of the complex takes place on crowded microtubules. For example (i) how does cargo-binding 
affect motor motility? and in connection with that - (ii) how does the presence of other motors 
(and also other cargos) on the microtubule affect the motility of the motor-cargo complex? 
We discuss these questions from a different perspective. The movement of a cargo is conceived 
here as a hopping process resulting from the transference of cargo between neighboring motors. 
In the light of this, we examine the conditions under which cargo might display bidirectional 
movement even if directed by motors of a single polarity. The global properties of the model 
in the long-time regime are obtained by mapping the dynamics of the collection of interacting 
motors and cargos into an asymmetric simple exclusion process (ASEP) which can be resolved 
using the matrix ansatz introduced by Derrida 



keywords - intracellular transport by molecular motors; bidirectional movement of cargo, 
traffic jam on microtubules; ASEP models. 

1 Introduction 

Research interest in the origins of the long-range bidirectional movement of particles (organelles, 
vesicles, nutrients) driven by molecular motors is motivated by fundamental questions concerning 
the nature of interactions between motors and their cargos as transport processes take place. A 
current explanation for the phenomenon relies on the idea that motors of different polarities act 
coordinately on the same particle at different times. If, however, they act in parallel, the bidirec- 
tional movement would reflect dominance of one or another kind of motor achieved by a tug-of-war 
mechanism [1], [2], [3], [5]. An important question that remains in this context concerns the 
mechanisms that would promote such coordination [6]. Alternatives to the coordination or tug-of- 
war models in the literature arise from the possibility of attributing the phenomenon to a dynamic 
role of the microtubules [7] or to a mechanical coupling between different motors [Sj. 

A general difficulty encountered within any of these views is related to the presence of other 
particles (including other motors) on the microtubule at a given time that are not directly involved 
with the transfer process. These other particles are expected to impose restrictions on motility 
and performance of the motors that are directly interacting with cargo at that time [9]. Contrarily 
to these expectations, however, data from observations of beads driven by kinesins in steady-state 
conditions indicate that the number of long length runs of such beads increases significantly as the 
density of motors at the microtubule increases, although their velocities remain essentially unaltered 
within a wide range of motor concentrations |T0], |Tl]. Thus, the reality of traffic jam in crowded 
microtubules still challenges the current view of long-range cargo transport that presupposes an 
effective and controllable movement of the motor (s) arranged into a motor-cargo complex. This, of 
course, requires a certain degree of stability of motor-cargo interactions and motor processivity. 

Our intention here is to discuss these problems from a different perspective by bringing into this 
scenario the model introduced in [12] to examine cargo transport as a hopping process. According 
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to that, motors and cargos would not assemble into complexes to put transport into effect. On the 
contrary, each motor would function as an active overpass for cargo to step over to a neighboring 
motor. In this case, the long-range movement of cargo is envisaged as a sequence of these elementary 
(short-range) steps either forwards or backwards. In [12] we examined the conditions under which 
this may happen, accounting for the fact that motor motility is affected by the interactions with 
other motors and with cargos on the microtubule. There, we considered the presence of a collection 
of interacting motors, all of them presenting the same polarity (kinesins may be thought of as 
prototypes) and a single cargo. Here, we examine whether it is possible to explain in a similar 
context the origin of the observed bidirectional movement displayed by cargos. 

The particular mechanism we propose to substantiate the hopping differs from that suggested 
in [12]. It keeps, however, the same general ideas of the original. As it will be explained below, we 
view the hopping of cargo between motors as an effect of thermal fluctuations undergone by motor 
tails. The flexibility of the tails may promote contact and, eventually, exchange of cargo between 
neighboring motors. 

As in [12j, the model dynamics is mapped into an asymmetric simple exclusion process (ASEP) 
[13] , [11] , [I5] whose stationary properties are resolved explicitly in the limit of very large systems. 
Other ASEP models have already been considered in the literature to study the conditions for 
motor jamming in the absence of cargo [9], [16], [I7]. Our model is conceived to account explicitly 
for changes in the dynamics of the motors that at a certain instant of time are interacting with 
cargos. 

The model is reviewed here in order to include a second cargo in the system, still keeping the 
presence of motors of a single polarity. We believe that this approaches more realistic situations 
in which the simultaneous presence of many cargos and motors on the same microtubule must be 
the prevailing situation [7]. We show that under these conditions, a cargo may be able to execute 
long-range bidirectional movement as it moves over clusters of motors assembled either at its back 
end or at the back end of the cargo in front. One may recognize in this a possibility for explaining 
the origins of self-regulation in intracellular transport since it has been suggested in the last few 
years that signaling pathways involved in intracellular traffic regulation can be performed simply 
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by the presence of cargos at the microtubule [18]. We then speculate that the passage of cargos 
on microtubules does not get blocked by motor jamming. On the contrary, jamming operates as 
an allied process to promote long runs of cargos across motor clusters. In this case, the density 
of motors on the microtubule can be identified as an element of control in intracellular transport 
since it directly affects the conditions for jamming. 

It is worth mentioning that the model developed here does not rule out other possibilities, 
such as the tug-of-war or competition models. What we suggest is that the presence of motors of 
different polarities may not be essential to explain the origin of the bidirectional movement. 

The hopping mechanism is presented in Sec. 2. The kinetic properties of the extended version 
are developed in Sec. 3, considering the presence of two cargos. In Sec. 4 we present our results. 
Additional remarks and conclusions are in Sec. 5. 

2 An alternative for cargo driving 

The stochastic model in [12] formulated in a lattice describes the dynamics of motors and cargos 
accounting for (i) steric interactions among different particles moving on the same microtubule; (ii) 
the presence of motors of a single polarity; (iii) the fact that cargos do not move if not driven by 
motors. 

The crucial point is in item (iii) because it requires a specific model for motor-cargo dynamics as 
transportation takes place. We offer here a slightly different view from that in [12] keeping however 
the reliance on the ability of motors to transfer cargo. One way by which this may be achieved is 
sketched in Figure (1). The stepping of cargo would be accomplished as it is released from a motor 
to which it is attached at a certain instant of time and then get attached to another (neighboring) 
motor either at the left or at the right - see Figure (la). This process, like the one discussed in |12] . 
relies strongly on the flexibility of the motor's tail. The idea was inspired by experimental results 
suggesting a dynamic role of the kinesin's coiled-coil segment in the process [19], [20] and also by 
data indicating that under load, kinesin motors display an oscillatory movement [21] [22]. Here, we 
think of these oscillations as signaling fluctuations in the position of motor's tail, not necessarily 
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being correlated with displacement of its center of mass. If this is the case, such oscillations would 
promote contact between neighboring motors favoring cargo exchange. 

Accordingly, long-range displacements of cargo would reflect a hopping process extended over 
many neighboring motors which may be accomplished if these motors get jammed into clusters for 
sufficient long periods of time. Notice that the whole mechanism does not require special stability 
of cargo-motor binding. On the contrary, the transfer of cargo by a motor would be ease by a loose 
attachment between them. 

The model in [12] was resolved explicitly considering the presence of a single cargo in the system. 
The averages for the quantities of interest were determined in steady-state conditions. We showed 
there that the long-range displacements of the cargo would occur predominantly in the backwards 
direction, i.e. in opposition to the direction of the movement of the considered motors. We shall 
show here that the same dynamics may lead cargo to display bidirectional movement if the system 
contains at least one more cargo interfering with the movement of the motors. 

To examine the properties of the system containing two cargos and an arbitrary number of 
motors, we map it into the same ASEP as in [12] whose dynamics can also reproduce the scheme 
in Figure ( 1 ). 

2.1 Movement of motors and cargos: the ASEP model with two cargos 

We consider a one-dimensional lattice with M sites, representing the microtubule with periodic 
boundary conditions. This system contains motors and a number K of other particles - the 
cargos - that interact with motors in order to move. Each site can be occupied by a motor or 
by a motor attached to a cargo (see Fig. 1(a)), otherwise it is empty. The total number of sites 
that remain unoccupied is G = M~N~K>{). Here, we analyze the long-time behavior of 
this system for K = 2 and determine the average cargo velocity as a function of the parameters. 
The results indicate conditions for cargos to perform a type of long-range movement that share the 
characteristics of the observed bidirectional movement. 

The map of the dynamics shown in Figure. 1 into the considered ASEP is carried out as follows. 
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First, each site is identified by its position j = 1, 2, ...N at the lattice. Then, to each of these sites is 
associated a variable aj that assumes integer values 0, 1 or 2 such that aj = if the site j is empty, 
aj = 1 if it is occupied by a motor; or aj = 2 if it is occupied by a motor attached to a cargo. 
With these, a configuration C of the lattice is specified by the set {ai(J2... ctn}- The dynamics of 
the ASEP that reproduces the elementary steps in Fig.l can now be defined. For this, consider 
that at each time interval dt a pair of consecutive sites, say j and j + 1 are selected at random. 
The occupancy of these two sites is then switched according to the following rules 

(a) 10 —7- 01 with rate k, probability kdt 

(b) 12 — )• 21 with rate w, probability wdt (1) 

(c) 21 — )■ 12 with rate p, probability pdt 

where the pair (j, j + 1) is represented by the values of the corresponding site variables {aj, crj+i).The 
parameters k, w and p are the assigned probabilities per unit time (rates) for occurrence of the 
processes indicated. Process (a) describes the possibility for a motor (kinesin) that carries no cargo 
to step forward to a neighboring empty site (Figure lb). Processes (6) and (c) account for the 
switching of the cargo between two neighboring motors. This accounts either for backward (6) or 
for forward steps (c). Notice that the dynamics conserves the number of particles of type-1 as well 
as those of type-2. 

In order to investigate the long-time dynamics of a cargo resulting from these elementary steps 
we use the matrix ansatz introduced by Derrida [23], [15]. The idea is to represent the probability 
Pn,m{C) of a configuration C of the system with sites and M particles of type-1 as a trace over 
a product of non-commuting matrices, each specifying the corresponding site occupancy: 

1 ^ 

Pnm^C) = Tr W{5„^^iD + 5,^,2^ + 5a.,o^) (2) 

' 1=1 

where 

N 

Zn,m = Y.'^r HiS^^^iD + 6,^,2A + 5.,,o^) (3) 

is the normalization. The sum runs over all configurations for which Ylf ^a-i,i = M and Ylf ^<Ti,2 = 
K = 2. In this product, a site i is represented by a matrix D if it is occupied by a motor (cij = 1) 
or by a matrix A if occupied by a motor with a cargo (cTj = 2); if the site is empty it is represented 
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by a matrix E (cTj = 0). In order to calculate averages over these configurations in the stationary 
state, it is necessary at first to find the algebra that must be satisfied by these matrices such that 
the probabilities defined in (|2]) satisfy the stationary conditions [Ti] . 



5^p^,M(c")r(c" ^ c) - p^,M(c)r(c ^c') = o (4) 

c 

where the sum extends over all configurations of M motors distributed over N — K lattice sites. 
Observe that the nonzero terms on the LHS of the above equation are those for which configurations 
C and C differ from each other at most by the positions of a pair of consecutive sites, which can 
be reversed by any of the elementary processes defined by the dynamics in ([T]). In this case, each 
factor r(C" —7- C) (or r(C — )■ C')) must be replaced by the rate w, k or p for the corresponding 
elementary process that brings C back from C (or C from C ). 

The algebra corresponding to the ASEP defined by the dynamics ([1]) has been presented in 
([12]) ioT K =1: 

DA - xAD = E-D 
DE = E 

(5) 

EA = E 
EE = E 

with 

X = (6) 

w 

Here, we shall use this same algebra to evaluate the traces over products of matrices D, A and E 
that appear in calculating averages over the quantities that characterize the movement of a cargo. 
Before proceeding, however, a few remarks are in order. 

2.2 Traffic profile in the system with two cargos 

The model with two cargos is not ergodic. The dynamics preserves the number of empty spaces 
in each of the two partitions defined by the initial positions of the two cargos in the system with 
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periodic boundary conditions. In this case, all configurations C and C'that satisfy equation (jl]) 
must share the number of empty spaces in each of the partitions. Moreover, configurations in which 
the empty spaces are all concentrated in one of the two partitions must be excluded, for these do 
not satisfy (jlj) with the algebra ([5]). 

We treat the initial conditions (/C), namely the number of empty spaces - h m one of the 
partitions and G—h in the other, as random variables. This artifact shall account for the uncertainty 
one has in experimental data regarding the relative positions of the particles, and also for effects 
of random processes that are not explicitly described by the present model such as motor binding 
and unbinding at the microtubule. For computing averages, we shall account first for all possible 
configurations at fixed h and then average the results over h. The procedure is further specified 
observing that (a) because there is no reason to favor any initial configuration, we may consider that 
h is uniformly distributed and (b) in analogy with a situation of equilibrium, we take the average 
annealing as the averages over particle configurations are performed in parallel with average over 
h. 

The measure PN,Mi.^(h)) of a configuration C^h) of a subset-h is written as the trace over a 
product of matrices A, D and E that satisfy the algebra in ([S]): 



In the expression above, each pi is a binary variable such that Pi G {0,1}, i = 1,2,...N and 
ki = 1 — Pi satisfying 

ki + k2 + ... + kk + kp+i + ... + k^-i + kN = h 



PN,M{C{h)) 



1 



(7) 



,{h) 

'N,M 



£)Pp JTI^p J^£)Pp+i ]TI^P+l £)Pn-i JT^f^N-i £)Pn JT^^n 



and 



(8) 



kk+2 + kk+3 + ■■■ + fcj 



'p ~ 



G-h 



The normalization 




(9) 
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is conveniently expressed in terms of the weier hts Wr ^ . ^. These are defined as 

the sum over the traces corresponding to the configurations that belong to the subset h for which 
the occupation of the sites ... i — + 1... and ...j — l,j,j + 1... in the n-tuples are fixed and 
specified by the values of the corresponding site variables (o".t_i, (jj, ...) and (crj_i, aj, <Jj+i, ■■■), 
respectively. 

The PN,M{C(h)) defined above must satisfy the stationary conditions 

PN,Mic[H))^ic[,^ ^ q.)) - PNMCih))ric^H) q,)) = o (lo) 

where the sum extends over all configurations C^'^^ that belong to the subset h. The transition 
rates r(C^x — )■ C(^h)) lead configurations C^n into configurations C(^hy 



2.3 The average velocity of a cargo 

Consistently with the above definitions we represent the average value < V(^h) > of the velocity 
of any of the two cargos at fixed h as 



< V{h) > 



7(h) 
'N,M 



(21)~(2) 



(12)~(2) 



The configurations associated with wl^2i)-(2) such that the specified cargo has one motor at its 
right side that allows it to move one step to the right at a rate p. Similarly, in the configurations 
associated to wl^^2)-[2) there is a motor at the left side of this cargo that allows it to move one step 
to the left at rate w. In both types of configurations the neighborhood of the other cargo is not 
specified. 

It shall be convenient to subdivide the above sum into sums over configurations having the same 
trace. This is achieved by specifying in f ITT]) the occupation of the sites that precede both cargos. 
For this, < V(^h) > is rewritten as 

1 



< 



V(h) 



> 



'N,M 



P 



lyih) _lU/W ^w'^''^ 

'^'^(121)-(12) ^ '^'^(121)-(02) ^ '^'^(021) -(12) ^ "^"^ (021) -(02) 



-W 



'^'^(12)-(12) ^ '^'^(12)-(02) 



(12) 



To proceed in the evaluation of (fT2|) it is also convenient to replace the site variables {cij} 
by block variables {m,} and {qi} i = l,2...k, that assume integer values to represent, respectively, 

sequences of motors and empty sites in a configuration C^h)- With these, the sum over configurations 

-(h) 

(121)-(12) 



that contribute to wj^l'L.^f,^. for example, can be expressed as 



{qi} {mi} 

with mfc+i, mfc+2, mp+i > 1. Here, D^^ (or E*^^) indicates a product of rrii (or qi) matrices D (or 
E). The symbol on the summation signals indicates that these are restricted to the configurations 
that satisfy the constraints in (|H]) for < h < G. 

All the traces in the RHS of f lT2|) can now be reduced with the aid of the algebra in ([5]). For 
this, we use the identity 

D^AE = x^AE (14) 
that also follows directly from [12j. The results are quoted as follows 

^a2)-(i2) =^'tr{ED'^^AED'^>=ED'^^AED'^^+^) ="^' x'^^x"'Hr{E) 

conf conf 

^(m)-(i2) = ^'tr(ED"^'y4D'^^-ED"^M) = ^'x"^^x'"^tr(E) 

conf conf 

^{S)-m =^'tr{ED"''AED"'^EA) =^'x"'Hr{E) 

conf conf 

<2i)-(i2) = E 'HEAD-^ED-^A) = J2 'x-Hr{E) (15) 

conf conf 

^a2i)-(02) = E 'tr{ED"^^AD"^^EA) = ^ 'x-^Hr{E) 

conf conf 

<2i)-(i2) = E 'tr{EAD-^ED-^A) = 'x-tr(E) 

conf conf 

<2)^(02) =Y.'t^{EAEA) =J2'tr{E) 

conf conf 

The above expressions are independent of h. The only dependence on h in the evaluation of 
the weights comes from the multiplicity of the configurations. Configurations for which h = or 
h = G contribute with a factor with respect to the contributions from all other configurations 
that result in the same trace. 
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2.4 Average over the random variable 



We now take the average of < V(^h) > over all realizations of h that assumes an integer value within 
the interval h E [1, G — 1] for G > 2, with equal probability. This is performed here as 

^ ^ EJP^(21)-(2) -t^H^(12)-(2)} ^ {pF(21)-(2) -wW(i2)-(2)} ^^^^ 

that corresponds to the average annealing in analogy to a situation of equilibrium. The averaged 
quantities are indicated by the bars over the corresponding symbols representing the weights W 
and normalization Z. 

Now, observe that because the traces do not depend on then the sums in the above expression, 
both in the numerator and in the denominator, account for all possible configurations of arbitrary 
sequences of empty and occupied sites, keeping fixed just the n-tuples indicated in each term. With 
this, the restrictions imposed on the sums in (ITSll are removed. 

2.5 Sum over configurations 

We estimate the number of configurations that contribute to W(^a-i.-^a-ia^+i )-(aj.i(Tj(Tj+i ) for a given 
n-tuple (crj_icrjcrj+i....) — (crj_ia"ja"j+i....) by fixing the relative position of the cargos ^ and counting 
for all possible sequences of O's and I's. We then sum over ^ observing the invariance of the trace 
under cyclic transformations. The results are compiled below. 
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^-^ ^-^ 2 \ M — rrij — rrii J 



0) W^(i2M02) - El M?^~ )[(iV-m,-3)x-]tr(E) 



mi=l 



e) W^(m)-(02) - E (!!_!'^._i)[(^-^' -4)x™»]Mi?) 



(/) W^(02iM02) ^ [(iV-4)]Mi?) 



(^7) W^(02)-(02) ^ 



M J 



Variables rrii and rrij indicate the number of possible consecutive motors at the left of the cargos 
in each of these configurations contributing to a given W[rTi.iaiai+i,,,,)-{aj^iajaj+i,,,,)- The sums over 
rrti and mj are estimated here in the limit of very large systems for which N ^ oo and M — > cxd 
keeping the motor density M/N — )■ p finite within the range < p < 1. In this limit, the sums 
converge to integrals and these integrals can be evaluated using Laplace's asymptotic method. 

Consider, for instance, the sum (a) in f[T7|) . We use Stirling's formula A^! ~ V^irNN^e'^ to 
approximate the factorials involving the variables N and M and define the new variables [23] 

z = rrii/N and y = rrij/N. (18) 
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z and y assume continuous values in this limit so that the referred sum converges to the integral 

4 fp l-(.p-y) /I _ „ _ r X^^X^y 



The function z) in the expression above depends only on the sum y + z : 



h{y, z) = h{y + z) = [l-iy + z)] ln[l - (y + z)] - [p - (y + z)] \n[p - (y + z)]. (20) 
Thus, by defining u = y + z, (fT9|) can be rewritten as 

where 

/2 = + z/lnx (22) 

In order to apply Laplace's method for estimating the above integral, it is convenient to change 
the order of the integration observing that 

rp rp rp 

/ dy / dui-) = du dyi-) (23) 

^0 Jy Jo Jo 

With this change, the integral in y becomes trivial and the double integral in fl2T!) reduces to 



/= fu. r^ ^ , e^^^(-)ci^ (24) 
Jo ^ P-^{l-yf 

which can be estimated by its maximum at large . For this, notice that f2{i^) has a maximum at 



_ 1 - xp 

1 — X 



Thus, if 



(A) xp < 1 (26) 

the condition z/max < P can not be satisfied and the maximum contribution to the integral in 
comes from the extremum of the interval at u = which is a local maximum. The result is l25l 



tA,C 

If ~ 

72 



I 1 exp[-iVplnp] ^2^^ 



2 J iV2 ^ [In(xp)] 
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If however, 



(B) xp > 1 (28) 
then z/max is locahzed inside the integration interval so that the integral in (HM is estimated as [25] 

{x - 1)2 1 (1 - xp) [2^ 



h2 ~ 



{xf {l-pf" 1 



X 




^exp{iv[lna;-(l-p)ln{fE3]}. (29) 



We use this same procedure to estimate all the remaining terms in expression f[T6|) . We merely 
quote the results below, making some extra comments when necessary. 

For estimating the sum indicated as (6) in (fTTll we notice that the difference Amj,mi between 
the expression in the RHS of (6) and the sum in (a) is 

\n„m, = Yl ^121-12 - Yl ^12-12 
conf conf 



mj,mi 



rtij ,mi 



[N - Mf ^'s^ (N - rrij -mi-?,\ x^^x^' 



' ^ \ M — rrij — rrii IN — m,- — m,- 

1 m,i=l ^ J « / 3 ^ 



(30) 



Applying Laplace's method to the resulting integrals after taking the thermodynamic limit, it gives 



m.j ,mi 



-7V(l-p)ln(l-p) 



2v/27riV(l -p)*' 

' 1 1 exp(-A^plnp) 
(Inxp) 



2 ' 



xp < 1 



or 



X < 



^fl 



\x 



5/2 



, (1 -a;)2 [2^ 



Inx - (1 -p) ln{fz^ 
xp > 1 



(31) 
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The sum over configurations that contribute to W^i2-o2 in (flTl) - (c) is estimated through the 
asymptotic behavior of a single integral, which gives 

N'^i^ — pY -Af(l-p) In(l-p) 



conf 



v/27riV(l - p) 

1 1 exp{—Np\np) 
N ^/p |lna;p| 

X < or 

1 {1-xf 



xp < 1 



exp < 



In X — (1 — p) In 



(1-p) 




27r 



(32) 



xp > 1 



Analogous procedures are used to estimate the sums over configurations of the kind W021-12 
and 1^121-02 in (fT7|) - (d) and (e), which coincide in this limit: 



-12 

conf 



1 exp(— A^plnp) 
|lnxp| 



-N{l-p) In(l-p) 



(33) 



X < 



or 



1 (1-a:)^ 

[ {xf'H^-pf" 



exp <^ N 



In a; — (1 — p) In -fe^ 



(£-1) 



xp < 1 



27r 




For the remaining sums indicated in ( fT71) - (/) and ((7), it is sufficient to estimate the relevant 
contributions to 0{-\fN) which are 

N {1-pY 



-02 

conf 



-N[{l~p) ln(l-p)+plnp] 



2 ^27riV(l -p)p 



(34) 



and 



-02 

conf 



AT P^'^P ^-jV[(l-p)ln(l-p)+plnp] 



(35) 



v/27riV(l-p)p 

In the following, we analyze the results for the average velocity of a cargo obtained in this limit 
using the estimates above. 
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3 Results 



The average velocity of a cargo in the system of interacting motors and cargos that obey the ASEP 
dynamics set in ([T]) can now be analyzed observing the differences in the expressions obtained above 
for the integrals in each of the asymptotic regions limited by the range of the product px of the 
two variables p and x. Such differences lead to distinct behaviors for < v > characterizing different 
phases of the system that, in turn, reflect the differences in the distribution of motors along the 
considered microtubule. 

A a fixed value of x such that x > 1 and for 



(a) < p < l/x (36) 
< f > is obtained from the behavior of the integrals for px < 1, resulting in 

{p-w) + \\Yipx\[2pp\\npx\ + {2pp-w)]-p{l- p) 

< V >~ 2 

1 + 4 |ln px\ + [In px] 
Within the complementary region in which 



(b) l/x < p < 1 (38) 

< f > is determined from the behavior of the integrals for px > 1; we find 

k 

<v> (39) 

X 

Notice that for any x < 1 the condition (a) px < 1 is always satisfied so that the results for 

< V > are given in this case by fl?r|) within the entire interval < p < 1. Thus, for small values of 
X the system does not exhibits phase transitions. 

The behavior of < f > is shown in Fig. 2 for the whole range of motor density p, at fixed k = 1 
and w = 3 and at various values of x ([6]). 

We notice in these results that for varying p and for x slightly above 1, the average velocity 
of the cargo changes sign. This means that at steady state, which may be achieved at sufficiently 
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short times after an eventual change in motor density at the microtubule, cargos may adjust and 
change their direction of propagation moving across motor clusters. 
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4 Discussions and additional remarks 



The mechanism for cargo transfer envisaged here is equivalent to a hopping process in which the 
associated rates depend on site occupation. Because motors move and their movement is affected 
by the presence of the cargos and all other motors on the microtubule, the long-time dynamics of 
the system must be examined globally. 

Our results indicate that the existence of mutual interactions and the fact that many cargos 
are allowed to coexist at the microtubule are determinant for reproducing in this context the 
characteristics of the bidirectional movement. We show that within a certain range of motor density 
a cargo in this system executes long-range displacements in both directions. We may argue then 
that long-range cargo transfer is facilitated by traffic and specifically, by the assembly of motors into 
clusters, which characterizes traffic jam. The presence of the other cargos in the system is essential 
for this to occur as they function as additional obstacles that interfere in the motor density profile. 
Each cargo induces aggregation of motors at its back end. In turn, this provides the conditions for 
cargo to execute long-range displacements either backwards, over the aggregate assembled at its 
back end , as well as forwards, over the aggregate assembled at the back end of the cargo in front. 

As it was originally formulated the model does not account for the possibility that a motor 
with one or more attached cargos may move as well. In fact, this is the only mechanism that 
is usually employed to describe cargo transport and it is the basis for coordination or tug-of-war 
models. Also, for simplicity we have considered interactions of a cargo with a single motor at a 
time. Should a set of motors be allowed to interact with the cargo to participate in the transfer 
process then the map into the ASEP would need to be modified accordingly. We are currently 
working on these possibilities by including into the dynamics ([T]) a process of the kind 20 — )■ 02 
that recovers ergodicity of the model [26] 

We should emphasize that the occurrence of the long-range bidirectional movement as a conse- 
quence of the hopping processes devised here may happen by the action of motors of just one kind 
possessing a well defined polarity. Changes in motor density and related traffic profile suffice as a 
mechanism to control cargo direction and the size of the runs determined essentially by the extent 
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of motor clusters at jamming conditions. This offers a rather straightforward explanation for the 
data mentioned above suggesting that the number of long run-lengths performed by the observed 
beads increases significantly as the density of motors at the microtubule increases [11]. In addition, 
the results presented here indicate that, for sufficiently high values of the motor density for which 
p > 1/a;, at the point where the model displays a phase transition, cargos would perform a uniform 
movement (on average) since their velocities become independent of p. Such behavior has also been 
observed in the same set of experiments. 

As noticed by Ma and Chisholm [6], "little is known regarding motor traffic and how it correlates 
with the movement of cargo" . Here, we offer a possibility based on the idea that the transport does 
not require the action of an external agent to coordinate the process, or a tug-of-war mechanism 
or even the existence of a mechanical coupling between two kinds of motors as proposed more 
recently |8]. Instead, it suggests that such coordination can be achieved by collective effects on the 
course of the dynamics as the system "self-organizes" so that it presents characteristics that reflect 
an internal (and global) order that does not have its origin in the characteristics of the external 
medium. It is then possible that the necessary transport in cells is accomphshed just by adjusting 
the density of motors at the microtubule. 

Accordingly, the presence of processive motors of different polarities that are normally required 
to explain the movement of a putative motor-cargo complex would not be necessary. Transportation 
here is based on a mechanism that requires formation of clusters of motors, not necessarily on their 
ability to travel along long distances. 

Acknowledgements 

I would like to thank Domingos H.U. Marchetti for very helpful and enthusiastic discussions 
regarding the conditions imposed by Fubini's theorem and the procedure used here to estimate the 
double integrals; and also Elisa T. Sena for pointing to me many of the difficulties in the initial 
developments of this work. I thank Scott Hines for kindly editing a preliminary version of the text. 



19 



This work had integral support from Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP) - Brazil. 



References 

[1] S. P. Gross, Hither and yon: a review of bidirectional microtubule-based transport, Phys. Biol., 
1, Rl-Rll (2004). 

[2] M. A. Welte, Bidirectional Transport Along Microtubules, Curr. Biol. 14, R525-R537 (2004). 

[3] K.B. Zeldovich, J-F, Joanny, J. Prost, Motor proteins transporting cargos, Eur. Phys. J. E 17, 
155-163 (2005). 

[4] R. Malhck, S. P. Gross, Molecular Motors: Strategies to Get Along, Curr. Biol., 23 R971-R982 
(2004); 

[5] S.Klumpp, M. J. Miiller, R. Lipowsky, Cooperative transport by small teams of molecular 
motors, Biophys. Rev. and Lett, 1, 353-361 (2006). 

[6] S. Ma, R. L. Chisholm, Cytoplasmatic dynein-associated structures move bidirectionally in 
vivo, J. CeU Sci. 115, 1453-1460 (2002). 

[7] I. M. Kuhc, A. E. X. Brown, H. Kim, C. Kural, B. Blehm, P. R. Selvin, P. C. Nelson, V. 1. 
Gelfand, The role of microtubule movement in bidirectional organelle transport, Proc. Natl. 
Acad. Sci. USA 105, 10011-10016 (2008). 

[8] S. Ally, A. G. Larson, K. Barlan, S. E. Rice, V. I. Gelfand, Opposite-polarity motors activate 
one another to trigger cargo transport in live cells, J. Cell Biol. 187, 1071 - 1082 (2009). 

[9] S.Klumpp, R. Lipowsky, Traffic of Molecular Motors Through Tube-Like Compartments, J. 
Stat. Phys. 113, 233 -268 (2003). 

[10] A. Seitz, T. Surrey, Processive movement of single kinesinss on crowded microtubules visualized 
using quantum dots, EMBO J. 25, 267 - 277 (2006). 

20 



[11] J. Beeg, S. Klumpp, R. Dimova, R. S. Gracia, E. Unger, R. Lipowsky, Transport of beads by 
several kinesin motors, Biophys. J 94, 538-541 (2008). 

[12] C. Goldman, E. T. Sena, The dynamics of cargo driven by molecular motors in the context of 
an asymmetric simple exclusion process, Physica A 388, 3455-3464 (2009). 

[13] P.A. Ferrari, C. Kipnis, E. Saada, Microscopic Structure of TraveUing Waves in the Asymmetric 
Simple Exclusion Process, Ann. Prob. 19, 226 - 244 (1991). 

[14] B. Derrida, S.A. Janowsky, J.L. Lebowitz, E.R. Speer, Exact solution of the totally asymmetric 
simple exclusion process: shock profiles, J. Stat. Phys. 73, 813 - 842 (1993). 

[15] R.A. Blythe, M.R. Evans, Nonequilibrium steady states of matrix-product form: a solver's 
guide, J. Phys. A: Math Theor. 40, R333-R441 (2007). 

[16] A. Parmcggiani, T.Franosch, E. Frcy, Totally Asymmetric Simple Exclusion Process with 
Langmuir Kinetics, Phys. Rev. E70 046101-1 - 046101-20 (2004). 

[17] Chowdhury, D. Traffic Flow of Interacting Self-Driven Particles: Rails and Trail, Vehicles and 
Vesicles. Physica Scripta, T106 (2003) 13 - 18. 

[18] K. J. Verhey, J. W. Hammond, Traffic control: regulation of kinesin motors, Nat. Rev. Mol. 
Cell Biol. 10, 765 - 777 (2009). 

[19] M. Imanishi, N. F. Endres, A. Gennerich, R. D. Vale, Autoinhibition regulates the motility of 
the C.elegans intraflagellar transport motor OSM-3, J. Cell Biol. 174, 931-937 (2006). 

[20] A.Yildiz, M. Tomishige, A. Gennerich, R.D. Vale, Intramolecular strain coordinates kinesin 
stepping behavior along microtubules. Cell 134, 1030-1041(2008). 

[21] M. Nishiyama, H. Higuchi, T. Yanagida, Chemomechanical coupling of the forward and back- 
ward steps of single kinesin molecules. Nat. Cell Biol. 4, 790-797 (2002). 

[22] Y. Ishii, Y. Taniguchi, M. Iwaki, T. Yanagida, thermal fluctuations biased for directional 
motion in molecular motors, BioSystems 93, 34-38 (2008). 



21 



[23] B. Derrida, M.R. Evans, The asymmetric exclusion model: exact results through a matrix ap- 
proach, Nonequilibrium statistical mechanics in one dimension. University Press, UK, Chapt. 
14, 277-304 (1997); B. Derrida, M.R. Evans, V. Hakim, V. Pasquier, An exact solution of a 
ID asymmetric exclusion model using a matrix formulation, J. Phys.A26, 1493-1517 (1993). 

[24] D.H.U. Marchetti, P.A.F. da Veiga, T.R.Hurd, The 1/N-expansion as a perturbation about 
the mean field theory: a one-dimensional fermion model. Comm. Math. Phys. 179, 623-646 
(1996). 

[25] J. D. Murray, Asymptotic Analysis, Springer (1984). 
[26] L. W. Rossi, T. Santos, C. Goldman, work in progress; 

[27] K. Kawaguchi, S. Uemura , S. Ishiwata, Equilibrium and transition between single- and double- 
headed binding of kinesin as revealed by single-molecule mechanics, Biophys. J. 84, 1103 - 1113 
(2004); 

[28] J. Howard, Mechanics of Motor Proteins and the Cytoskeleton, Sinauer Associates Inc., 2001. 



22 



Figure Caption 



Figure 1 - Dynamics of motors and cargos. (a) Cargo transfer. It liappens here through 
a mechanism of hopping between neighbor motors. Due to the fiexibihty of the tail, the attached 
cargo may display small oscillations leading to the possibility of it being caught either by the 
motor at its left or by the motor at its right. The corresponding processes 12 — > 21 or 21 — )■ 12 
are represented in the figure, (b) The step of a motor. The time spent by the motor with the 
two heads attached to the microtubule is much larger than the time it spends with just one of the 
heads attached ^7], as a part of the "hand-over-hand" mechanism proposed to explain the kinetics 
of two-headed motor proteins [28J. Occupation of a site by a motor occurs here whenever it is 
occupied by the two heads of a motor. The motor step is then represented as 10 — 01 which is 
indicated in the figure. 

Figure 2 - The average velocity of a cargo for k = 1 and w = 3 as a function of the density p 
of motors at the microtubule, for various values of parameter x. 
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